Spin-Hall effect in triplet chiral superconductors and graphene 
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We study spin-Hall effects in time-reversal symmetry (TRS) broken systems such as triplet chiral 
superconductors and TRS preserved ones such as graphene. For chiral triplet superconductors, we 
show that the edge states carry a quantized spin-Hall current in response to an applied Zeeman 
magnetic field B along the d vectoa-, whereas the edge spin-current for B _L d is screened by 
the condensate. We also derive the bulk spin-Hall current for chiral triplet superconductors for 
arbitrary relative orientation of B and d and discuss its relation with the edge spin-current. For 
TRS invariant system graphene, we show that the bulk effective action, unlike its TRS broken 
counterparts, does not support a SU(2) Hopf term but allows a crossed Hopf term in the presence 
of an external electromagnetic field, which yields a quantized bulk spin-Hall current in response to 
an electric field. We also present an analytical solution of the edge problem for armchair edges of 
graphene and contrast the properties of these edge states with their time reversal symmetry broken 
counterparts in chiral superconductors. We propose possible experiments to test our results. 

PACS numbers: 73.43.-f, 73.61.Wp, 74.20.Rp, 74.70.Pq 



I. INTRODUCTION 



The generation and control of spin currents in various 
solid state systems have been an area of active research 
interest. Such systems can be broadly grouped into two 
categories; those which have broken time-reversal sym- 
metry (TRS) and those which do not. Examples of the 
first group include chiral superfluids and superconduc- 
tors such as two dimensional (2D) ■^Hc-A films^, 2D lay- 
ered strontium ruthenate^ and d + id superconductor"^. 
In these systems, the superfluid or superconducting 
order parameters break TRS and also provide a gap 
for the quasiparticle states at Fermi surface. Such a 
gap, together with broken TRS, naturally lead to chi- 
ral edge states and Hall effect^*^. However, none of 
these systems support the usual charge Hall effect, since 
they are either charge neutral ("^He-A) or do not have 
conserved physical charge density (ruthenates, d + id 
superconductors)-"'^. Instead, they support spin and 
thermal Hall currents which were found to be quan- 
tized with N = 2 in d + id systems^, and A^ = 1 in 
^He-A and triplet superconductors^. The second class 
of systems which do not break TRS but exhibit spin- 
Hall effect include doped semiconductors with spin-orbit 
interaction's and grapheneii. In graphene, which pro- 
vides a realization of Haldane's modeli^, the spin-orbit 
interaction provides the gap for the effective low-energy 
Dirac-like Hamiltonian of the electrons near the Fermi 
{K and K') pointsii. Such a gap breaks TRS individ- 
ually for each spin species, but leaves the whole system 
invariant. Consequently, one expects a quantized charge 
Hall current js in the presence of an external electric 
field separately for the spin species s =ti-l- These cur- 
rents move in opposite directions so that j-\ — ~ji . Thus 
the total charge Hall current vanishes; however, there ex- 
ists a net spin-Hall current jspin = U] ~ Ji)^/2e. In the 



absence of a Rashba term in the Hamiltonian, this spin 
current is exactly quantized with a spin-Hall conductiv- 
ity as — e/27r. In the presence of the Rashba term, where 
the number of up and down spin electrons are not con- 
served, the spin-Hall conductivity deviates from its quan- 
tized valueiS. This bulk picture of the spin-quantum Hall 
effect in graphene is further supported by numerical veri- 
fication of the presence of chiral edge states at the bound- 
aries of the 2D graphene sheet Uii^. It has also been 
shown that although the spin-Hall conductivity deviates 
from its quantized value, the spin-Hall effect and the edge 
states are robust against weak Rashba couplingiii^. 

The quantization of spin-Hall conductivity in both 
classes of such quasi 2D systems necessitates the exis- 
tence of topological terms in their low energy bulk effec- 
tive action. The nature of these topological terms de- 
pends crucially on whether the system is invariant under 
TRS or not. For systems with broken TRS, it has been 
shown that the low energy effective action has a SU(2) 
Hopf term which may lead to quantized spin-Hall current 
in the presence of an external magnetic field gradient^iSiS. 
However, the details of the spin response of the triplet 
superconductors is expected to be qualitatively differ- 
ent from both singlet d + id superconductors where the 
pair potential does not break spin-rotational symmetry, 
and "^He-A, which, being charge neutral, do not expe- 
rience Meissner screening. A detailed understanding of 
the spin-current in triplet superconductors therefore re- 
quires an analysis of its effective action in the presence 
of an arbitrary Zeeman field. Such a study has not been 
undertaken so far. In contrast, for the second class of 
systems which respects TRS, one expects a crossed Hopf 
terniiS. which leads to a spin-Hall current in the presence 
of an external electric field. However, such a term has 
not been explicitly derived for graphene. 

The presence/absence of TRS is also crucial for the 



properties of the chiral edge states in these systems. Such 
states are well-known to exist for both triplet and sin- 
glet {d + id) chiral superconductors^i^i^. In these sys- 
tems, the edge modes corresponds to localized Bogoli- 
ubov quasiparticles states which have zero charge, but 
carry a finite charge current^'^. For singlet d + id su- 
perconductors, it is well known that the edge states also 
have a definite spin quantum number and would there- 
fore carry a quantized spin-Hall current in the presence 
of an external magnetic field^i. However, the spin struc- 
ture of the edge states in triplet chiral superconductors 
has not been studied before. In contrast, for systems 
which preserves TRS, the edge modes, in the absence of 
external perturbations, can not carry any net charge cur- 
rent, but is expected to carry a finite non-quantized spin 
currenliiiii^. However, as we shall show, they can carry a 
experimentally detectable quantized charge current, if a 
population imbalance of spin-up and spin-down electrons 
is created at the edge. 

The main results reported in this work are the follow- 
ing. First, for triplet chiral superconductors with the 
pair-potential given by 



A(k) = Ao (s • d) (fcj, - iky) /kp 



(1) 



where s denotes Pauli matrices in spin-space, kp denotes 
the Fermi wavevector, and d denotes a direction in spin 
space orthogonal to the spin of the Cooper pair«i, we 
show that the spin structure of the chiral edge states is 
qualitatively different from the previously studied d + id 
superconductors^. In particular, these edge states carry 
a quantized spin-current in response to an external mag- 
netic field B only if the external magnetic field is ap- 
plied along the d vector. The spin-current for B _L d 
is screened by Meissner current of the condensate. Sec- 
ond, we obtain an effective action for the bulk system 
in the presence of an arbitrary external Zeeman field, 
and derive an expression for bulk spin-current from this 
effective action for arbitrary relative orientation of the 
applied Zeeman field B and d . We compare this bulk 
spin-current with its edge counterpart for both B || d and 
B _L d. In the limit of zero external Zeeman field our ef- 
fective action reduces to those derived in Ref. 8 9 and 
contains an SU(2) Hopf term. Third, for the TRS invari- 
ant system graphene, starting from the low energy Dirac- 
like Hamiltonianii, we derive an effective action in the 
presence of an external electromagnetic field, and show 
that it contains a crossed Hopf termi^ which leads to a 
quantized spin-Hall current in response to an external 
electric field. We also include a weak Rashba term in the 
HamiltonianiiiiSi and demonstrate that the spin-current 
deviates from its quantized value in the presence of such a 
term. Finally, we present a solution of the edge problem 
in graphene for the armchair edge and obtain analytical 
expressions for the energy dispersion and wavefucntions 
of the edge states. The properties of the edge states ob- 
tained from this analytical derivation matches previous 
numerical resultsiiii^. We also point out that these edge 
states can carry a quantized charge current in response 



to an applied magnetic field which can be experimentally 
measured. 

The organization of the paper is as follows. We dis- 
cuss the properties of triplet chiral superconductors in 
Sec. m In Sec. Ill Al we develop the edge state picture 
and discuss the spin structure of the edge states. This 
is followed by Sec. IIIBI where we obtain a bulk effec- 
tive action for triplet chiral superconductors and derive 
an expression for the spin-Hall current in these systems. 
We compare the bulk spin-current with their edge coun- 
terparts for both B II d and B _L d. Next, in Sec. IIIII we 
study spin-Hall effect in graphene. In Sec. IIII Al we de- 
rive a bulk effective action for graphene and show that a 
crossed Hopf term exists in its effective action which leads 
to an quantized spin-Hall conductivity in the absence of 
Rashba coupling. We also compute the deviation of this 
conductivity from its quantized value in the presence of 
a weak Rashba term in the Hamiltonian of the system. 
Next, in Sec. IIII 51 we present an analytical solution for 
the edge state spectrum for the armchair edge and discuss 
the properties of these states. This is followed by a dis- 
cussion of possible experiments and concluding remarks 
inSec.HVI 



II. CHIRAL TRIPLET SUPERCONDUCTORS 

For triplet superconductors. Cooper-pairing occurs in 
the L — I, S — 1 channel, and consequently the pair- 
potential (Eq. ^ breaks spin-rotational invariance. This 
is manifested in the presence of the d vector in the ex- 
pression of the pair-potential (Eq. ^ which refers to a 
direction orthogonal to the direction of the spin S of the 
Cooper pair, so that a choice of spin quantization axis 
along d implies opposite-spin pairing^. In the presence 
of such a pair potential, the action for chiral supercon- 
ductors can be written as 

S = f d^rdti:1{t,r) idt - Mei^iX/) - Ep) 

-r+A(k)-r_A*(k)]^(t,r), (2) 

where r denotes Pauli matrices acting in particle-hole 
space, e(k) = (fc^ -I- ky)/2m is the kinetic energy, V = 
{ips, gss''4'^ " ), for spin indices s, s' —'\i, is the four com- 



ponent Pauli spinor, and g^ 



i{sy)ss' is the metric 



tensor in spin spaced. Notice that under time-reversal 
symmetry A{kx,ky) -^ A{—kx,ky) and thus the pair- 
potential breaks time reversal symmetry. We shall use 
this action to analyze the edge and the bulk properties 
of the system in Sees. Ill A\ and IHBI 



A. Edge states 

In this section, we analyze the spin properties of the 
edge states. We shall assume that the d vector is fixed 



along a specific direction due to the spin-orbit interac- 
tion throughout the superconductors. This is expected 
to be the case for layered chiral superconductors Sr2Ru04 
where the d vector is fixed along c-axis perpendicular the 
Ru — O planc^®. In this section, we shall adapt the fol- 
lowing convention for the choice of spin-quantization axis 
for the electrons: if an external magnetic field is applied, 
the spin-quantization will be chosen along the field; else 
it will be chosen to be along d. For the rest of this work, 
we shall choose h = c = 1, unless explicitly mentioned. 

First, we consider the edge problem in absence of the 
Zeeman field. Here we shall consider a semi-infinite sam- 
ple with an edge at a; = and assume a step-function 
dependence of the pair-potentialii. In this case, the 
Bogoliubov-de Gennes (BdG) equation can be written 
as 



Ho A{kF)s, 
A*(kF)s, -Ho 



*(a;,/ci.y) = E'^ix^kpy) 



(3) 



where Tig = ^i—idx, ky) — Ep = —ivpdx is the linearized 
dispersion and kpy is the in-plane momentum component 
of the Fermi surface along the edge, and ^P (x, kpy) = 

I Us {x, kpy) , i{(Ty)ss'v'' {x, kpy) ) Is the four-component 
spinor wavefunction. It is well-known that in the pres- 
ence of an edge which enforces a boundary condition 
^ (x = 0, kpy) — 0, Eq. 13 admits localized chiral edge 

statea^i^i^^i^'^ 



^ {x,kpy) — 2i^/K,e 



ikpyVp — KX 



sin {kpxx) 



sgn(s)wj 



E{kpy) = 



Aok 



Fy 



K — - Kpx, 

kpVp 



(4) 



where sgn(s) = ± for s =][, Us = exp{i9i,)vs, Ob 
is the global U(l) phase of the pair-potential A, and 
Ufi are the coefficients of spin-up and spin-down com- 
ponents of the normalized single-particle wavefunctions 
with I Ms I + \vs\ = 1. These chiral edge quasiparticles 
are BCS quasiparticles with equal proportion of spin s 
electrons and spin s holes; consequently, they do not 
carry any charge but have a well defined a spin quan- 
tum number sgn(s)/2 along d. They have a group ve- 
locity Ve = dE (kpy) /dkpy = Ao/fci? and carry a charge 
currenliSiiSiiS 
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(5) 



where Us (kpy) — 6 [kp — kpy) is the density of states for 
edge state with spin sgn(s)/2 at zero temperature. The 
charge current is not carried with group velocity and is 
not quantized. The presence of such a charge current 
at the edge of a superconductor also leads to Meissner 
screening current by the condensate, so that the sum of 
magnetic fields due to both the currents vanish at the 
bulk of the superconductoriSiiSiiS. 
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FIG. 1: Chiral edge states in spin triplet superconductors in 
the presence of a Zeeman magnetic field B along d. The red 
(blue) solid lines denote occupied spin-up (spin-down) edge 
states while the corresponding dotted lines show the unoccu- 
pied states. The dashed black lines are guide to the eye. 



In the absence of any external Zeeman magnetic field, 
the spin-up and spin-down edge quasiparticle states are 
equally populated and hence the spin-current carried by 
them cancels. Now, let us apply a Zeeman magnetic field 
B along d so that (choosing the spin-quantization axis 
along B and d) the BdG equations for the edge quasi- 
particles states can be written as 



A*(kF)s. 



= E'i'{x,kpy)l6) 



where /is is the Bohr magneton and g denotes the gy- 
romagnetic ratio. We note that here the effect of the 
magnetic field is to shift the energies of the spin-up 
and spin-down states in opposite direction: Eg (kpy) = 
Aokpy/kp — sgn{s)gfj,BB. As a result, there is an im- 
balance in the population of the spin-up and spin-down 
quasiparticle states as shown in Fig. ^ which leads to a 
net quantized edge spin-current 



(TZ\spin 

V-'ll^dgc 



1 



Ve{N^~N^) 



1 



-gUsB. 



(7) 



Note that the spin current, in contrast to the charge cur- 
rent, is associated with the group velocity v^, of the edge 
states. Further with our choice of spin-quantization axis 
(along d) , the condensate comprises of Cooper pairs with 
opposite spin pairing. Hence the condensate do not carry 
any spin or spin-current and can not screen the edge spin- 
current. 

Such a quantized spin current in the presence of ex- 
ternal Zeeman field is also carried by edge states in sin- 
glet d + id superconductors regardless of the direction of 
the applied magnetic field. The situation for the triplet 
superconductors is different as we now show by studying 



the spin response of the edge states in the presence of the 
Zeeman field apphed perpendicular to d. Here, without 
loss of generality, we choose d = y and take the spin- 
quantization axis along H = Hz. The BdG equation can 
now be written as 



"^0 - gfJ'sBsz 



* (x, kf 



ACkp) Sy 

= E'9{x,kFy).{8) 



One immediately notes that in this case, the mag- 
netic field does not shift the energy of the edge states, 
but shifts the Fermi wave-vector for the spin-up and 
spin-down quasiparticles in opposite directions: kps = 
\JkF + sgn{s)2mg fiB B . The BCS quasiparticles are now 
equal superposition of electron and hole states with the 
same spin s. As a result, just as in the case of charge, 
they do not carry a spin quantum number along z: 

= 0. The edge states now carry 



\ charge 



a net charge current 
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/edge 



E(^^ 



-^charge 



kFs 



hk 



kj='Tj—(} 



m 



"■Fs 



Snm 



(9) 



spin, the Meissner current generated as the response to 
edge charge current also carry a net spin-current 



xspm 
Meissner 
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(12) 



Thus the edge spin current generated by the response of a 
magnetic external Zeeman field B _L d is screened by the 
condensate, analogous to the charge current. Therefore 
the response of the edge states to an applied Zeeman field 
in triplet chiral superconductors depends crucially on the 
direction of the applied field. We shall discuss possible 
experiments to probe this behavior in Sec. IIVI 

Before ending this section, we note that all of the re- 
sults derived can be applied to ruthenates where the pair- 
potential is expected to have horizontal line of nodes with 
the substitution A ^- A |cos(cA:^)p". In particular, it can 
be shown that the spin response of the edge states remain 
qualitatively similar in the presence of such horizontal 
line of nodes. 



Now let us consider Meissner screening of this charge cur- 
rent. Since we have chosen the direction of spin quanti- 
zation axis perpendicular to the d vector, we have equal 
spin pairing in the bulk (A|| = A^^ = A) and hence 
the Cooper pairs carry a net Sz — ±?i. In the pres- 
ence of a Zeeman magnetic field, the spin-up and the 
spin-down condensates see a shifted chemical potential 
/i|l = M + ±ff/^B^- In this situation, one can apply the 
arguments of Refs.lsBlSin a straightforward manner and 
show that the edge currents would be screened by spon- 
taneously generated Meissner current 



A^ 



r 



= -{IsA 



\ charge 
/edge ' 



(10) 



so as to cancel the net magnetic field due to the edge 
currents in the bulk. 

The edge states now also carry a net spin current, since 
the electron and holes with spin s constituting the BCS 
quasiparticles have opposite velocities. In the presence 
of an external Zeeman field B _L d, kp-^ and kpi are 
different, and consequently the net spin-current is given 
by 



/jz ^spin 



edge 



-E' 



?n(s)(/. 



\ charge 
/edge 



9fJ.BB 

Stt 



(11) 



Notice that the spin current is not carried by the group 
velocity. Further, since the condensate now carry a net 



B. Bulk effective action 

To derive the bulk effective action, we begin with the 
action of chiral triplet superconductors in the presence of 
an external Zeeman field J5(i,r) and an arbitrary d field 
configuration. In the presence of such a term the action 
can be written as 

-T3(e(-z9, - Cr*) - Ef) - T+A(-a, - Cr') 
-r^A*{-zd,-Cr)U{t,r), (13) 



where Cq^^ = gfiB^ ■ B(i, r) is the Zeeman field and A is 
the pair-potential defined in Eq.^ Here, and in the rest 
of this section, we shall use the notation 9^ = {dt, dx, dy) 
and adapt the sign convention that all covariant vec- 
tors have F^ — {Ftj—Fxj—Fy) and all operators have 
dfj_ — {dt,dx,dy). The corresponding contravariant vec- 
tors and operators can be obtained by applying the met- 
ric tensor g^^, = g^^" — (1,-1,-1). Also, here we have 
followed Ref. |a to treat the Zeeman magnetic field as 
the time component of an external SU(2) gauge field 
(^ext ^ -^/2s„ • rj;j°'^* with Q.f^ = g^iBB. The spatial 
component of CJj'^* are fictitious fields which shall be set 
to zero at the end of the calculation. The advantage of 
introducing these fictitious gauge fields becomes appar- 
ent when we note that the spin-current can be obtained 



using these fields as 



SScs 



Sflf 



(14) 



where S(,g is the effective action obtained by integrating 
out Fermions in 5* (Eg. lT^ . 

Next we introduce a SU(2) rotation in spin space 
through the transformation ilj{t,r) -^ i/j'{t,r) = 
W^]{t, r), so that the local SU(2) rotation matrix rotates 
the d vector to z: U'^ {s ■ d)U — Sz- One can now write 
the action as 



S = Sq + Si 



So 
Si 



f27r' 



3 ^'T(fc)Go V/(fc), 



(15) 
(16) 






x^{afc,Go-i(fc),C^°*-'(p)}^V'(fc), (17) 



where the notation {..} , denotes anticommutation, the 
Green function GQ{k) is given by 

Go(fc;Ao) = - + -3(e(k)-i..) + Ao(k.r).. ^^^^^ 



uj'^ - El+ m; 



E^ = J {e{k) - EfY + Al 



and the gauge fields (7*°*^' are given by 



(19) 



c 



'total 






y^int 



C" 



-zc/ta^[/ = is„f7«'"*. 



(20) 



Note that here the internal gauge fields C' are 'pure 
gauge fields in the sense that they satisfy 



d^n 



int 



dM 



int 



rj!"* X n 



int 



/r-o. (21) 



These fields are small when the d fields are slowly vary- 
ing. Thus for a slowly- varying configuration of both d 
field and external Zeeman field B , one can carry out a 
gradient expansion in C'°*'*'. To this end, we now inte- 
grate out the fermions fields from 5 and carry out a gra- 
dient expansion of the resultant effective action in (7*°**^^ 
as done in Refs. |3|9 22. We shall concentrate here only 
on the terms in the effective action which contributes to 
the spin-Hall current. These diagrams are shown in Fig. 
121 After some straightforward, but tedious algebra, one 
obtains the effective action which has the same form as 




FIG. 2; Feynman diagrams which contribute to the SU(2) 
Hopf terms effective action. The closed circles represent 
1/2 {OGq^ /dki, Si} with b = fj,,!/ or X and i — a, /3 or -y 
as appropriate. The space-time indices /i, v and A take values 
t,x and y whereas the spin indices a, (3 and 7 tafce values x, y 
and z. The straight lines denote the Green function Go and 
the wavy lines denote the total SU(2) gauge fields ^2'°'"' 



that obtained for ^ He-A in Ref. 



'5'cff 


— 


*5int + "Jcxt + "-"cou 


pling, 






(22) 


'S'int 


= 


^e^^x 1 d^dtni'^^'dxnf^'''', 




(23) 


•^oxt 


= 


,^,^.^^ J d'rdt 


'V ■'■''A 










-in-*.(f2rxf^r)' 


7 




(24) 


'^coupling 


^ 


,11 ''^'^ 


Nef,^x 










x{d,d-ni-'xd)-{r,fxd) 










+2P {dtd - fij:''* X d) • (f^;;' X 


d) 


.(25) 


e f "''* = 


d^fll^' - 9,f2^'^' 


-ft 


_=xt ^ f^oxt ig 


the field- 



strength corresponding to the external SU(2) fields, e^i,A 
is the antisymmetric tensor and the coefficients N and P 
can be written in terms of the Green function Go^k; Ag) 

8 Q 



N = 2ttTt 



Go(fc;Ao) ^^°"lf'^°^ Go(fc;Ao) 



dky 



dkx 

-Go(fc;Aoj 



,aGo-\fc,Ao)^ ,,.^ ,9G-i(fc,Ao) 



dbJ 



P = 27rTr 
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dkr 



Go(fc;Ao) ^^°''f'^"^ Go(fc;Ao) 



Go{k;-Ao 



dGo\k,0) 



dky 



(26) 



(27) 



where Tr denotes matrix traces and sum over all frequen- 
cies and momenta. Notice that N is precisely the topo- 
logical invariant that characterizes the SU(2) Hopf term 
as shown in Refs. I8l9l For chiral triplet p-wave supercon- 
ductors and •^He-A, it is well-known that TV = 1. The 



coefBcient P, however, is not a topological invariant and 
its numerical value depends on the details of the system. 
Evaluating the integrals and matrix traces in Eq. [23 one 
finds that P = 1 as long as we restrict ourselves to clean 
systems and T = 0. In this work, we shall restrict our- 
selves within this domain. We also note that N is exactly 
quantized only if there are no nodes in the superconduct- 
ing gap; however, the deviation of N from its quantized 
value is expected to be vanishingly small in ruthenates in 
the presence of the line of nodes, since the gap vanishes 
only along a line of the entire three dimensional Fermi 
surface. We shall ignore such deviation in the rest of this 
section. 

In the absence of external Zeeman field, we are left with 
only the internal fields JTJf^In this case, the effective 
action S^s = Sint (Eqs. I22I23|) reduces to SU(2) Hopf 
term derived in Refs. laS The contribution to the spin- 
current, on the other hand, comes from S'coupUng and S'cxt, 
and thus requires the full effective action in the presence 
of a finite Zeeman field. Using Eqs. El |^ and 1251 we 
obtain the spin-Hall current for arbitrary configuration of 
the d(r, t) and for any relative orientation of the applied 
Zeeman field B and d 



where (J'lpcd'JJ; i^ the edge spin current given by Eq. 
H Thus the bulk spin current is equal and opposite 
to the edge current. This has an interesting conse- 
quence. Imagine that we place the semi-infinite sam- 
ple with in a solenoid coil centered around the edge of 
the sample at a; = which produces a magnetic field 
along d. In this case, the solenoid will create a both a 
bulk and the edge spin current so that total spin cur- 
rent {q)Z"c + (^y)buik = 0. This is due to the fact 
that the solenoid creates a potential difference but no 
net electrochemical potential difference in the sample. 
Consequently, any net transverse spin-current must van- 
ish. This effect is analogous to the well-known absence 
of charge quantum Hall current when a semiconducting 
sample is placed between the plates of a capacitor. Thus 
we conclude that, contrary to the claim in the Ref. ^, the 
bulk spin current does not vanish for B || d; however the 
bulk and the edge spin current do cancel out each other, 
and the total spin-current vanishes. Also note that in 
this case the orbital effect of the magnetic field, which 
we have neglected so far, do not contribute to the spin- 
current since the condensate do not carry any net spin. 
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(28) 



In contrast, when the magnetic field is applied perpen- 
dicular to the d vector (we choose B = Bz and d || y), 
the bulk spin-current becomes 



where i and j takes values x,y, 7 = g^s, and Cy is the 
antisymmetric tensor with exy — 1. Note that Eq. [5H| 
has to be supplemented with the equation governing the 
dynamics of the d fields in the presence of an external 
magnetic field. This has been derived, in the presence of 
a pinning term, in Rcf. 2^ For the rest of this section, we 
shall concentrate on the case where the pinning potential 
(and hence the pinning frequency) is larger than all other 
scales in the problem, so that d is fixed along the c-axis. 
First, let us consider the case when B |j d = Bz. In 
this case one obtains a quantized spin-Hall bulk current 
from Eq. Ei 

II N 1 

on OTT 

This result differs from the conclusion of Ref. la, where 
the bulk spin current is claimed to vanish for magnetic 
field along the d vector. To understand the reason for 
this difference, let us consider a superconducting sample 
occupying semi-infinite space x > with a edge at x = 0. 
In this case, since magnetic field B decays to zero deep 
inside the superconductor, the bulk spin-current is given 

by 



(I, 



y /bulk 



1 

'8^ 



OvYdx 

{gf,BB) = -m)Zl^ (30) 



i^y )bulk 



(Jyfdx 


P ^ X 1 , 

— {gfiBB) = —{gfiBB). 



(31) 



Thus we find that there is a net bulk spin-Hall current 
for B _L d which is in accordance with the conclusion of 
Ref. 8. Notice that in this case, as we have seen in the 
Sec. IIIAI the edge spin current is already screened by 
the Meissner current and hence it does not contribute to 
the net spin current. By taking into account the genera- 
tion of Meissner current, we have qualitatively taken into 
account the orbital effect of the applied magnetic field. 
A detailed numerical computation, following the lines of 
Ref. 5, is required to address this issue on a more quanti- 
tative level and is beyond the scope of the present work. 
Also, we note that (/y)buik depends on P and not the 
SU(2) Hopf coefficient N, as can be easily seen from Eq. 
I2H1 Thus we conclude that the bulk spin-Hall conductiv- 
ity in a triplet chiral superconductor (and also ^He-A) for 
B _L d is not truly quantized, although for clean systems 
at T = where P = 1, it reaches its quantized value 
I/Stt. 



III. GRAPHENE 

The low energy Hamiltonian of graphene, in the pres- 
ence of Rashba couphng, is given hy^ 

H = j d?ril:\v) - ivp {r^a^d^ + <7ydy) + AiT^a^^s^ 

(32) 





where vp is the Fermi velocity, r^, ct^ and Si are the 
Pauli matrices in the K — K', A — B and spin spaces 
respectively, Ai is the spin-orbit coupling gap, and Xr 
denotes the Rashba coupling. Note that the spin-orbit 
coupling gap does not break time reversal symmetry on 
the whole, but does so for each individual spin species. 
Here "0 is a 8 component fcrmionic spinor field given by 

V- = (v-f T , V'f T , V'f i , V'f j , V'Ai > ^Bi > i'Ai > V'f i) • In the 
rest of this section, we shall analyze this Hamiltonian in 
the presence of an weak external electromagnetic field, 
and also study its energy spectrum in the presence of an 
armchair edge. Our derivation can in principle to other 
models for spin-Hall effect which has similar structure of 
the effective low-energy HamiltonianSi 



A. Bulk effective action 



In this section, we derive a low energy effective action 
of graphene in the presence of an external electromag- 
netic field and obtain the expression of spin-current from 
that effective action. First, we study the system in the 
absence of Rashba coupling. We begin with the action 



S = j Srdti^Hr.r) ^^, - eA, - v,{r.aA-^^. - eA.) 



ay{-idy - eAy)) - AiTza^s-d{r,t) 



^(r,T),(33) 



where A^ denotes the [/(I) gauge fields corresponding to 
an external electric field Ei = —dtAi — diA^. Note that 
here we have introduced an unit vector field d(r, r). The 
Hamiltonian H ('Ea. l32|) corresponds to the configuration 
d = z and we shall set d = z at the end of the calculation. 



Next, as in Sec. IIIBI we introduce a SU(2) rotation in 
the spin-space V'(r,T) — *■ ^'(r,r) = U'^iI){v,t) where U is 
a local SU(2) rotation matrix which rotates the local d 
vector to z: U^ s ■ dU = Sz- The partition function can 
now be written in terms of the new spinor fields ip'ik) = 



FIG. 3: Feynman diagrams which contribute to the possible 
Hopf terms effective action. The left panel contains diagram 
contributing to S^^' and S^^' where as the right panel shows 
diagrams contributing to S" . The filled squares represent 
the vertices dG^^ /dkt with b = fi or v, as appropriate. The 
closed circles represent 1/2 {9G^^/9fct,, s^} as in Fig. ??, 
whereas the straight lines denote the Green function Go. 



Z = f ViP'^V^p'VCf.VAf.expiiiSo + Si)] , 
So = I ^'^i''Hk)Go\k)^'ik), (34) 

^1 = - / '^■zrr'" i^'Hk+p) 



d kduj 
"(2^ 

/d^kd^pdwdpQ ,,4., 

x^{9/c,Goi(fc),e^4p) + C^(p)}^V'(fc),(35) 

where the fields C^ — —iU'^d^U = 1/2 s • il^ are the 
SU(2) gauge fields, the configuration d = z corresponds 
to Cp — and we have used the same notations as in 



Sec. Ill Bl Here the Green function Ga{k; Ai) is given by 

uj + vpiTskxaoo + kyCFy) + AiT^a^Sz 



Go(fc;Ai) 



El 



IT] 



Ek = JvUkl + kD + Al 



(36) 



Note that the tp' fields see a gap term AiGzTzSz and 
therefore describes the same action as the one corre- 
sponding to H (Eq. I32|l in the absence of the C^ fields. 
The advantage of introducing the d fields now become 
clear, since we find that the spin-current density j" = 

V''^^ {9fe^G(^^(fc), s"} V') can be obtained as 



dSas [Cn , A., 



af7f 



(37) 



nf=o 



where S'es is the effective action obtained by integrating 
out the Fermion fields. Here the /z = t component of the 
spin-current is simply the spin-density j( = (ip''^ sip'') . 

To obtain the effective action, we now integrate out the 
fermionic fields and compute the effective action for A^ 



and C^ fields. Since for slowly varying d field configura- 
tion C^ is small, we can carry out a gradient expansion 
in powers of the fields A^ and C^ and their derivatives in 
a straightforward manner. The relevant terms in the ef- 
fective action which contributes to the Hopf terms comes 
from the first derivative of the polarization bubble and 
triangle diagram shown in Fig. |21 These terms are given 

by 



coefficients can be written as 



5(2) 



^(3) ^ 



5«+5(2)^^(3)^ (3g) 

a^,x fd'rdtA^dxA^, (39) 

+Cl j d'rdtn'^.^^^nl (40) 

d^^^^j d\dtA,dxn^^. (41) 



where we have used the definition C^ = \/2soS^'^, and 
all repeated indices are summed over. The coefficients 
of the terms in the effective action can be expressed in 
terms of Green functions as 






xGo(fc;Ai)— ^Go(fc;Ai; 
dkx 



(42) 



C^>({^-^}/o(.;Ar){^, 



S/3 



xGo(fc;Ai)— i^Go(fc;Ai; 
dkx 



+ 
(43) 



jdG^ 



-1 



xGo(fc;Ai)— ^Go(fe;Ai- 
dkx 



(44) 



(45) 



where Tr denotes trace over all Pauli matrices and inte- 
gration over frequencies and momenta. We now need to 
evaluate these coefficients and we begin with h°'^^ and 



h^P 



We note that for b J^y^, all the coefficients with 
a ^ z and (3 — z ot vice versa vanish when we take 
trace over Pauli matrices in spin-space. The same line of 
reasoning shows that only the coefficients c""^ or c^^^ 
can be non-zero. Further it is easy to see that only d^ ^^ 
has non- vanishing trace. Thus the only non- vanishing 



-1-1 



dG~^ 



^-(-^4Trr^C„(.;-A.)^Co(.;AO 



^^^"(''^^))=<"^' 



(46) 



dG. 



5Gn 



Tr ^.^n!^Go(fc; _Ai)^j^Go(fc; A,) 



^"^ 2"^ V"^ dkt, """''"' ""' dk^ 

x^^Go(fc;Ai) 

dkx 
xxz _ yyz _ _uyx 



(47) 






dG~^ \ 

x-^Go{k;Ai)j ==a^^A/e^ 



(48) 



rf^.A = -^TrL^Go(fc;Ai)^Go(fc;Ai) 



x^G„(fc;A,)) = 



dk 



'-'fj.i^X- 



(49) 



a/3 



It turns out*- that coefficients b J^-^ also satisfy the iden- 
tity b"l^ — -^b"^-^. Using this, we find that the contribu- 
tion of terms proportional to b^J^^^^ in the effective action 
leads to a total derivative which we ignore and those due 
to c""^ vanish identically. Further one can use the fact 
that the fields J7^ are pure gauge fields and hence satisfy 
the identity 

df,n^ - d^^Qf, -n^xn^ = o. (so) 

Using Eqs. 1201 and 07| it is easy to show that all the 
terms proportional to b^Jlx vanish. Thus finally we are 
left with the effective action 

S^^^ = a^,xjd\dtA,dxA^, (51) 

I. (52) 

5(3) = dl^xjd\dtA,dx%. (53) 

The next task is to evaluate the coefficient a^^Xi which 
can be done by straightforward evaluation of integrals. 
But before resorting to algebraic manipulation, it is eas- 
ier to note that one needs a^,yx = e^iniyxa for electromag- 
netic gauge invariance of S^^' , where a is given by 



^(2) ^ ^ d\dtnidxn 



dGn^ 



5^-(l7°"<'^^'>^°""^^'' 



x^Go(fc;Ai) 



(54) 



Further, we find that the contribution to a in Eq. [SH 
is a sum of contributions from spin | and | sectors: 



a — a-[+ai. Whereas each of these two contributions can 
be finite since the Hamiltonian fEa l32|l breaks time rever- 
sal symmetry for each spin species, their sum must vanish 
since the total Hamiltonian for both the spin species is 
time reversal invariant. Thus we conclude that a — Q and 
hence 6^^_j^ = 0, so that graphene does not support charge 
Hall effect in the absence of an external magnetic field 
and also does not have a pure SU(2) Hopf term. However, 
the same argument tells us that d^^^^ = —ee^^^xd ^ 0, 
since d — a-^ — a-^^ 7^ 0. Indeed a straightforward evalua- 
tion of a shows flj = — a| = — l/27r. Hence we conclude 
that the effective action of graphene supports a crossed 
Hopf term 



9(3) - —f ^ 



d^rdtVLlF^ 



(55) 



where F^x = duA\ — d\Ai, is the electromagnetic field 
tensor. This leads to a quantized Hall spin-current (Eq. 

EZIi 



Jt 



'2tt 



-E. 



'3' 



(56) 



where E — —WAq — dfA is the applied electric field. 

Next we introduce the Rashba term in Eq. |221 We 
shall restrict the analysis to the case of A/j/Ai ^ 1 and 
assume that the strength of the Rashba term is not suf- 
ficient to destroy the spin-Hall phase. In the presence of 
such a term, the inverse of the Green function in Sq (Eq. 
EJ becomes 



g; 



-IR 



-Afl {r^Syax - SxO-y) . (57) 



One can now repeat the same analysis as described 
above with Go(fc; Ai) replaced by G^{k\ Ai, \r). Since 



9Gn ^ /dk^ is identical to 9Gq ''/dk^, one finds that all 

l^ap a/37 



the arguments regarding the coefficients a^^A, ^u^a' '^u^> 



and d" A remain the same and, in the end, we are left 
with 



5(3) 
d^ 



d"e 



I^LUX 



d^dtniF, 



uXi 



2 r^^fc^^°^^'^^'^^^^fe 



y 
xG^ik; Ai, X^)^!f±—Gll{k; A,, Xj,) , 



B. Edge states 

The edge states in graphene in the absence of any ex- 
ternal magnetic field can be analytically obtained from 
the Dirac Hamiltonian H fEa. l32|l for the armchair edge. 
Throughout this section, we shall restrict ourselves to the 
case Afl = 0. The geometry we study here is that of a 
semi-infinite sample occupying x > with an armchair 
edge at a; = 0. For such an edge, the boundary condi- 
tion demands both ipA and iJjb to vanish at a; = 0. In 
what follows, we construct an analytical solution for the 
low-energy subgap localized energy states of the graphene 
hamiltonian H fEa. l32l) which respects the boundary con- 
dition ^Pa{x — 0) = and ^Pb{x = 0) = 0. Notice that 
in doing so by starting from Dirac equation, we ignore 
the lattice effects at the edge, which usually puts a mo- 
mentum cut-off kc above which the present description 
breaks down. The estimate of kc can not be reliably 
done without going into numerical analysis of the lattice 
Hamiltonianii*24. 

To construct such a solution, we note that the only 
way to achieve the above-mentioned boundary condition 
is to superpose wave-functions in the K and K' points. 
We therefore try a normalized wavefunction of the form 



V^edge y-^j '^y ) 



'""v^ 



,K 



,K 



m- 



.K' 



,K' 



iK'-: 



(59) 

Here s =1, | denotes the spin of the electrons, K and K' 
are the wavevectors of the K and K' points in Brillouin 
zone of graphene, k~^ denotes localization length of the 
edge states, and ky denotes momenta parallel to the edge. 
The normalization condition '^ g J dx \ipcdgcix)\' 



the wave functions imply '^^ u^^ 

and the boundary condition '0cdgc(a^) = 
condition 



K{K') 



= 1 of 

2 

= 1 



''As 






necessitates the 



(60) 



From Eqs. 021 and EHl it is easy to see that 
the wavefunctions {u^st'^Bs) ^^d (w^g,u||^j satisfy 
Schrodinger equations 



27r 



6A2 



(58) 



where the ellipsis represent terms which are higher order 
in A/j/Ai, and we have obtained the last line by exphcitly 
evaluating the matrix traces and frequency and momen- 
tum integrals in the expression of d^. Thus we find that 
the spin-Hall conductivity deviates from its quantized 
value and equals e(l — X%/{6AI))/2t: for Xr/Ai <C 1. 



sgn(s)Ai ivp{K — ky 
ivpin + ky) — sgn(s)Ai 



-sgn(s)Ai —ivpin + ky) 



-ivpin 



sgn(s)Ai 



'■As 



'■Bs 



E 



,K' 



,K' 



'As 
'Bs 



= E 



,K' 



,K' 



(61) 



where sgn(s) = +(^) for s —] {[). Solving the 
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FIG. 4: Edge states near an armchair edge in graphene. The 
sohd and dotted hnes show occupied and empty spin-up (red) 
and spin-down (blue) edge states. An equal occupation of 
the spin-up and spin-down edge states result in cancelation 
of charge current. 



The edge states, in contrast to their chiral supercon- 
ducting counterparts, have opposite group velocities for 
spin-up and spin-down electrons as shown in Fig.^ Fur- 
ther, in a transverse momentum state ky, the edge states 
carry a charge current 

jyih) = evFn{ky) {ipcdgd o-y lipcdgc) ^ esg-a{s)vFn{ky), 

(68) 

where we have neglected the small 0{k/\Kx — K'^\) con- 
tributions. Note that the charge current is carried with 
the group-velocity Vg = sgn{s)vF and thus has opposite 
direction for opposite spins. Hence these edge states do 
not carry a net charge current for equal occupation of 
spin-up and spin-down electron states, but carry a net 
spin-current circulating along the edge. Note that an 
estimation of the total spin-current carried by the edge 
requires a knowledge of the states at large ky , and hence 
can not be reliably done within this formalism. How- 
ever, we would like to stress that the low-energy linear 
dispersion of the edge states with opposite group veloci- 
ties for spin-up and spin-down electrons agrees well with 
the numerical calculations of Refs Ji^^. 



Schrodinger's equations, one obtains 



IV. DISCUSSION 



^K E-sgn{s)jE^ +vI{k'^ -k^) 



,K 



'■Bs 



ivpin — ky) 



E + sgn(s) ^E^ + vl{K^-kl) 



ivpin + ky) 



VFK ^ J/^l-E^ + vlkl 



(62) 

(63) 

(64) 



The boundary condition ■(/'(x = 0) = (Eq. I60|l now 
yields the condition 



£;-sgn(s),/£;2-t-4(K2_fc2) 



E + sgn{s)jE-^ + vl{K-^~kl) 



Ky ~\ ri 



(65) 



Together with Eq. 1641 this gives us localized states at 
the edge with energy E{ky) and localization length k^^ 



E{ky) = sgn{s)vFky, k ^ 



VF 



(66) 



V'edge(a;,fcy) = ^^^l-yV-^^e-^-^r.(s)./i 



zsgn(s) 



-,JKr 



isgn(s) 



JVi'r 



(67) 



Notice that the analysis of the Dirac Hamiltonian pre- 
dicts linear dispersion of edge states for all ky. This 
is clearly an artifact throwing out the lattice and so 
we would expect the results to break down at momenta 



^max = /j^^ ag commented earlier. 



The most direct way of detecting spin/charge current 
at the edge is to create a population imbalance of the 
spin- up and spin-down edge states. Such experiments 
has been performed for quantum Hall samples in Ref. 
I25I and suggested for ruthenates in Ref. 18. The idea is 
to apply a magnetic pulse B(t) to the edge of a chiral 
triplet superconductor with B || d. This creates a imbal- 
ance in the population of the spin- up and spin-down edge 
states leading to local magnetization which travels with 
the group velocity Ve- The magnetic moment due to this 
spin imbalance can be detected by a squid magnetome- 
ter as pointed out in Ref. 18. The spin structure of the 
edge states can also be verified since a similar experiment 
with B _L d will produce a null result. This feature dis- 
tinguishes the triplet chiral superconductors from their 
singlet counterparts where the direction of the applied 
pulsed magnetic field does not matter. 

A similar experiment can be performed to verify the 
edge state picture for graphene. Let us consider creat- 
ing a population imbalance for sub-gap edge states in 
graphene by applying a magnetic field pulse B along z. 
In the absence of spin- flip scattering and for S < Ai 
where the states in the bulk are gapped, such an im- 
balance would persist and lead to a circulating quan- 
tized charge current. The spin- flip scattering rate in 
graphene is expected to be small since the Rashba cou- 
pling Afl ~ 0.5mK is small compared to the spin-orbit 
gap Ai ~ 1.2K, and thus we expect the population im- 
balance and hence the charge current to persist for a 
sufficiently long timeS&. The magnitude of this charge 
current depends only on the population imbalance be- 
tween up and down spin low-energy subgap electrons and 
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hence, unlike the edge spin-current, can be rehably esti- 
mated from the edge state picture developed in Sec. 1111 1j1 
Using Eq. EHl we obtain for the magnetic pulse along z 



7-chargc 
edge 



evF'^{n^{ky) ~ ni{ky)) 



2e 



5Ms^-(69) 



This quantized charge-current for gi-J-sB = Ai = 1.2K 
leads to a net current of I^^^l^'^ — 8nA which can be 
easily detected by a setup similar to that proposed in 
Ref. |23- One of the central feature of the edge states in 
graphene is that the direction of this charge edge cur- 
rent depends on the direction of the applied magnetic 
pulse. This feature can also be detected by a time-of 
flight measurement analogous to those suggested in Refs. 
Il§l25l Alternatively, one can also selectively populate the 
up or down-spin edge states by spin-polarized tunneling 
into the edge states. The direction of the edge current 
measured in such a tunneling experiment would depend 
on the polarization of the tunneling electrons and would 
serve as a definitive proof of the spin-dependent velocity 
of the edge states. 

In conclusion, we have studied the bulk and edge prop- 
erties of TRS broken and TRS invariant systems such as 
triplet chiral superconductors and graphene. We have 



shown that both the structure of the edge states and 
the topological terms in the bulk effective action of these 
systems depend crucially on whether TRS is preserved 
or not. For TRS broken systems such as triplet chiral 
superconductor, the edge states carry a net charge cur- 
rent and the bulk action contains a SU(2) Hopf leading 
to a spin-Hall current in response to the gradient of ap- 
plied Zeeman magnetic field. We have also shown that 
both the edge and the bulk spin-Hall current of triplet 
superconductors, in response to an applied Zeeman field 
B, depend crucially on whether B || d or B _L d. This 
is a consequence of broken spin-rotational symmetry in 
these systems and serves as a distinguishing feature of 
triplet chiral superconductors from their singlet counter- 
parts. In contrast, for TRS invariant systems such as 
graphene, the edge states carry a net spin current while 
the bulk action contains a crossed Hopf term leading to a 
spin-Hall current in response to an external electric field. 
We have also suggested experiments to verify our some 
of our results. 
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